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The cranked Skyrme III effective Hamiltonian is applied for the analysis of the rotational depen- 
dence of the quadrupole and octupole moments in Ra, Th, and U isotopes. A comparison of the 
intrinsic electric dipole moments calculated in the model with available experimental and theoretical 
values is presented. It is found that the non-axial octupole deformation Y32 becomes favorable at 
. . . high spins for the actinide nuclei. 
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. In the last decade extensive experimental data have been obtained on low-lying states of negative parity (see for 

Q-C review ||l|,||). For the Ra - Th{Z - 88, iV - 134) and Ba - Sm{Z - 56, TV - 88) nuclei low 3" states, parity 

, j^ , doublets, alternating parity bands with enhanced dipole transitions have been found. Since the octupole interactions 

are strongest when pairs of orbitals from the intruder sub-shell (l,j) and the normal parity sub-shell (I — 3, j — 3) are 

00 I near the Fermi surface, the octupole correlations are expected to become important for these nuclei ^^■ 

^^ ■ The features observed in nuclei are very similar to ones familiar from molecular physics. In molecules a stable 

octupole deformation leads to the appearance of the rotational bands with the alternating parity levels connected 

by strong El intra-band transitions. MoUer and Nilsson used the assumption about a symmetry breaking of the 

ly-s ■ intrinsic reflection symmetry of the single-particle potential and demonstrated the instability with respect to the 

QQ , axial octupole deformation in the Ra-Th region Isl . The analysis was performed within the macroscopic- microscopic 

^^ ■ approach that combines the liquid drop model and the Strutinsky shell correction method for a modified harmonic 

0^ ' oscillator potential. Later, within this approach numerous results, based on various single-particle phcnomclogical 

^^ • potentials, have been reported regarding the presence of a stable axial octupole deformation in the ground state of 

some nuclei in the actinide region pi . The macroscopic- microscopic approach was also used in the combination with 

the cranking model to analyze properties of rotating nuclei with a stable axial octupole deformation (see e.g. ||6|,|7|). 

The problem of non-axial octupole deformations, being for a long time a terra incognita, attracts more attention 

during last years. The study of the non-axial octupole deformations could shed light on the tendency for the system 

^ ' on the way toward fission to avoid the superdcformation H. Calculations by the macroscopic- microscopic method 

^ with the Woods-Saxon potential predicted the importance of the banana-type Y31 deformation for highly deformed 

^ nuclei |8|, while the 132 deformation was found to be important in the ^^^Ra nucleus |H|. Manifestations of pronounced 

'^ shell effects have been discovered when non-axial octupole deformations are added to a harmonic oscillator model 
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Although the macroscopic-microscopic approach is a very efficient practical tool for the analysis of shell effects 
in finite Fermi systems, basically, it is a phenomenological method which is unable, for example, to predict binding 
energies. Development of microscopic approaches, which start with the effective nucleon-nucleon interactions treated 
within a self-consistent Hartree-Fock (HF) method, opened a new avenue in the nuclear structure studies (0,0. Using 
the Skyrme III effective interaction, Bonche et al. |ll| did the first self-consistent calculations for the ground states 
of the Ra-Th region. It was found that the total energy has a minimum for the nonzero value of the Q30 momentum 
in ^"^"^Ra. Robledo and Egido [Q employed the Hartree-Fock- Bogoliubov method with the Gogny interaction and 
found the axial octupole deformation in the ground states in ^^^•^^'^Ra, "^"^^Th and ^'*^~^^^i3a. Recent self-consistent 
HF+Skyrme calculations [Of suggest that the oblate states in nuclei with A '-^ 80 are soft against the Yaa deformation 
in the ground states. The self-consistent cranking calculations based upon the HF-I-Skyrme |1J] and the HF-I-Gogny 
p5[ approaches predict that the Y^i deformation is important for a correct description of the yrast band of ^'^S. 
However, the systematic investigation of the non-axial octupole deformations in the actinide region at high spins is 
missing. On the other hand, the calculations in the framework of the spdf IBM model show that the 1 —pf boson limit 
(dipole-octupole vibrations) can describe essential features related to the octupole bands in the Rn-Th region p6[ . 
Using the quasiparticle spectrum generated by the mean field with the reflection symmetry, the quasiparticle-phonon 
model reproduced quite well the experimental data for i?(i?l)-transition probabilities in ^^^"^'^^Ra and 229-23irp]_^ 
|l7| . As is seen from this discussion, the question of the existence of stable octupole deformations in the actinide 
region at different angular momenta is still controversial and challenging for microscopic approaches. In this paper 
the dipole, quadrupole and octupole intrinsic moments of rotating nuclei from the Ra-Th-U region are studied within 



the HF+Skyrme approach. To this aim we exploit the novel version (1.75) of the code HFODD |1^ for the Skyrme III 
force [Q in the particle-hole channel. The code was successfully tested for the analysis of superdeformed nuclei (see, 
for example, |2^). Therefore, the basis employed by the code is sufficient for our study. In the present calculations 
the pairing correlations are neglected. They are not expected to play an important role at large rotational frequencies, 
however, they may influence results obtained for the ground state. 

We briefly remind the main ingredients of the calculations performed with the code HFODD, since it will be 
useful for our discussion below. The code |21| solves the Hartree-Fock problem in rotating frame with the cranking 
Hamiltonian in the form 

Hn=T + V + V,oui -niy- J2 ^-^- + Y. Cam(Qam - Qa^)'- (1) 

T=n,p A/j 

Here T represents the kinetic energy operator, V is the Skyrme potential, Vcoui stands for the Coulomb term. The 
quantity f2 is the rotational frequency, ly is the total angular momentum projection onto the cranking axis (y- 
axis); A^. and Nr are the chemical potential and the particle number operator, respectively, for protons (t — p) and 
neutrons (r = n). The last term in (m) represents multipole constraints. It was introduced in order to find a good 
minimum of the total energy that corresponds to the expected symmetry of the mean field (see below) . Qxp, is the 
mass-multipole- moment operator (with the multipolarity A and the projection /z), Q\^ is the expectation value of 
the corresponding multipole moment for a given rotational frequency f2. The stiffness parameter C\^ determines 
the weight (or our priority) of the constraint for a given A and /i during the minimization process. The detailed 
description of each term in (|l|) can be found in |lqj2l[| . Minimizing the expectation value of the cranking Hamiltonian 
Hn, < HF; il \ Hq \ HF; fl > , for a given rotational frequency O and a fixed set of the Skyrme parameters, the yrast 
state and equilibrium values of mass-multipolcs momenta, < HF; fi | Qxf^ \ HF; Q > can be found. The knowledge 
of possible symmetries of a searched state could simplify the minimization procedure. We used the following options 
of the code: 

i. The mean field is supposed to be symmetric with respect to the transformations Sy = Pe^"^«, Sj — rPe^'"^=", 
S'J = TPe"^'^^' (T is the time reversal operator, P is the parity operator). Such a situation corresponds to the 
standard i?-invariance H when the expectation value, < HF; 17 | Qa^ | HF; fi >, is nonzero only for even values 
of A and /i. These symmetries are characteristic of well-deformed nuclei without any octupolc deformations. 
There are three symmetry planes (y — z,x — z,x — y) of the mean field. 

ii. The only mean field symmetry is Sy = Pe^™^y (there is only one x — z symmetry plane). This is a general case 
when all expectation values < HF; il \ Q\f^ \ HF; fl > can be nonzero. 

In the case i) the quantum numbers signature r = ±i (an eigenvalue of the operator Ry ~ e^*'^^») and parity tt = ±1 
can be ascribed to each HF single-particle state. In the case ii) the parity is not a good quantum number and the 
HF single-particle state can be characterized by the quantum number simplex: s = ±i (an eigenvalue of the operator 
Sy) (see a discussion about different symmetries in |22|]). In each cases discussed above one can reduce the number 
of iterations in the HF problem by the appropriate constraint value Qa^ m the last term in (n^). In addition, the last 
term in (|l|) allows to avoid physically meaningless "deformations" which could correspond, for example, to a shift of 
the nucleus center of mass. In order to ensure the translation invariance of the Hamiltonian we used the constraints 
with Qio = Qii — and the constraint Q21 = Q2-1 = for the principal-axis condition. 

The calculations for three chains of isotopes (Ra, Th, and U) were performed without imposing a restriction on 
the axial reflection symmetry (the option ii)). For each nucleus two minima of the total energy were obtained: one 
corresponds to nonzero octupole moments and the second minimum is related to the mean field solution with the 
conserved reflection symmetry (the option i)). In order to obtain the reflection asymmetric minimum the optimal 
values for the constraint Qso = 10'^ fm^ and for the stiffness parameter C30 ~ 0.01 were found and used in the 
calculations. The constraints for the other octupole moments result in the smaller binding energies. Both the minima 
of the total energy are very close. For instance, for a well octupole deformed nucleus ^^^Ra the octupole deformed 
minimum is lower by 1 MeV with respect to the reflection-symmetric minimum. All plots presented on Figs. 1-4 are 
related to the reflection asymmetric minima. From these figures it follows that the octupole deformed Ra, Th, and U 
isotopes have axial symmetric shapes at small rotational frequencies fl. With the increase of the rotational frequency 
both the non-axial moments Q22 and Q32 increase whilst the axial moments Q20 and Q^q decrease for all three elements 
studied. The exceptions are provided by the lightest isotopes of Ra and Th (^^^~^^"i?a, ^^^~^^°T/i). In these isotopes 
the moment Q30 reaches the maximum at n '^ O.SAIeV and then decreases. Similar tendency regarding the Qso 
moment has been predicted for even-even rotating isotopes 220-228^^],^ ^^ ^^iq niacroscopic-microscopic approach [pSJ. 
However, in the calculations p3] only axially symmetric shapes of the Woods-Saxon potential have been considered. At 



larger rotational frequencies, we obtained that heavier isotopes, especially, of uranium, behave as non-axial octupole 
deformed systems. In Fig. 3 the intervals of the rotational frequencies il for which the reflection asymmetric minimum 
lies below the reflection symmetric one are denoted by arrows. From Fig. 4 it follows that the non-axial octupole 
moment Q32 increases with the increase of the rotational frequency in ^^^'^^^'^^^Ra, 223,229,23irp]^ ^^^^^ 226,231,233-238^ 
isotopes. This predominance of the non-axial octupole deformation in some nuclei may be tested by measurements of 
the El transitions. A small displacement between the center-of-mass and the center of charge, caused by the non-axial 
deformation, may enlarge the El transitions within a parity doublet. In Fig. 5 the induced electric dipole moment 
Do, 

N Z 

Do = e[^Qw{p) ~ jQioin)], (2) 

is shown for different mass numbers A and different values of the rotational frequency. In (0) the quantities Qio{p) 
and Qioin) are the proton and neutron contributions to the mass-dipole-moment Qio, respectively. As is seen in Figs. 
3 and 5, the value of Do is correlated with the axial octupole moment Q30. The larger the value of Q30 the larger 
is the electric dipole moment Dq induced by the octupole deformation. This tendency is consistent with all known 
results concerning the Do (see M and the references there). The non-axial octupole deformation Q32 influences the 
value of Do as well, however, its contribution is less essential. In the chain of isotopes, from the lighter to the heavier 
one, we obtained the sign change of the quantity Dq^ which preserves at different rotational frequencies. This behavior 
is similar to the one predicted in |Q with the Gogny forces for Ra and Th isotopes without rotation. A small value 
of Dq preserves to high rotational frequency in ^^^^^^^Ra, 226-228rp]^ g^^^^ 232,233y_ Small changes in the quantity Dq 
over the full frequency range is obtained in ^^^~^^'*Ra, 226-228^^]^ ^^^^ 230-233u^ Notice that for ^^^"^^"'Ra our results 



are consistent with the experimental observation of the behavior of Dq |24|. The noticeable changes in the Dq due 
to the rotation are obtained in 228,23i,233,234j^j^^ 220,230,23i,233,234Th and 5^-227u. in Fig. 6 the calculated values of 
Do I are compared with the results of the other approaches and with the experimental values extracted from the 
i3(£'l)-transition probabilities M. We remind that these experimental data provide only the absolute value | i^o I- 
The comparison is done for two values of the rotational frequency ft = OMeV (low spin region) and fi = Q.lMeV 
(high spin region). One can see from Fig. 6 that our results are in agreement with the other calculations and that 
they reflect a qualitative behavior of the experimental data for Ra and Th isotopes. It would be interesting to test our 
prediction for U isotopes within other approaches and by experimental measurements. In our analysis the quantity Do 
is nonzero and has a physical meaning only in the intervals denoted by the arrows in Fig. 3. However, the correlations 
caused by residual interactions can also contribute to the values of the total energy p^ and multipole moments [p6| . 
Consequently, although the main features should be preserved, the range of the intervals may be modified. The work 
in this direction is in the progress. 

In summary, we analyzed the mass-multipole-moments Q2p, {fJ- — 0, 2), Qs^ (/i — 0, 1, 2, 3) and the intrinsic electric 
dipole moment Do in the Ra-Th-U region. The analysis has been done within the self-consistent cranking+HF 
approach using the Skyrme III effective interaction. It was found that the octupole moment Q30 decreases whilst the 
moment Q32 increases with the increase of the rotational frequency Q. A similar behavior is found for the quadrupole 
moments Q20 and Q22- The quantity Do changes the sign at il = OMeV a.t A ^ 224 (for Ra), A ^ 227 (for Th), and 
at A ~ 232 (for U). This tendency preserves over the whole range of values of the rotational frequency 0.0 < fl < 0.6 
MeV. 
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Figure Capture 

Fig.l The quadrupole moment Q20 as a function of the rotational frequency. The results are displayed: for Ra- 
isotopes at the top panel, for Th-isotopes at the middle panel and for U-isotopes at the bottom panel. 

Fig. 2 The quadrupole moment Q22- Similar to the Fig.l 

Fig. 3 The octupole moment Q30. Similar to the Fig.l. See also the text. 

Fig. 4 The octupole moment Q32. Similar to the Fig.l 

Fig. 5 The electric dipole moment Dq induced by the octupole deformation as a function of the rotational frequency. 
Similar to the Fig.l 

Fig. 6 Experimental (with error bars) and calculated intrinsic dipole moment | Dq \. Experimental and calculated 
values are indicated by symbols denoted: "els" and "tls" for measurements and our results for low-spin states 
(/ ~ Oh), respectively; "ehs" and "ths" for measurements and our results for higher spin states (/ ~ 8?i), 
respectively. Experimental values are taken from |0|. The symbols "BNls" and "BNhs" correspond to the 
results from [E3| obtained at fi = OMeV and fl = O.lMeV, respectively. The calculated values denoted by 
"ER" are taken from fill. 
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